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Errata 

Pages 10 and 11 : The captions for table 3 and t&ble 3 concluded 

should read: 

dY 

THE FUNCTIONS - £ FOR AIR ( Y = 1.4) FOR SEVERAL VALUES 

Pk dT 

OF THE INDEX k 

Pages 12 and 13 : The captions far table 4 and table 4 concluded 

should read: 

jy 

THE FUNCTIONS - =£ FOR AIR (y = 1.4) FOR SEVERAL VALUES 

pk dr 

OF THE INDEX k 


Addenda 


For completeness the additional formulas of reference 1 used 
in' the computation of the functions are as follows: 

For arbitrary positive indices 

I] C ( T ) = ^ k+lj rj 

Far negative nonintegral indices 

\(r) = t* F(a ]f *, Vk, 1-k; t) 

where 

a k + ^k = k ” P 

= - 1 (k + !)P 


and 


w _ \ >1 ab a(a 

F(a, b, c; t) = 1 + — r + — 


l)b(b + 1 ) 
2 ! c(c + l) 
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Far negative integral Indices 

(a k - k) (-b^ - k) (a k - k) (a t - k + l) (n k - k) (* k - k + l) 

^ ' " 1 ! (k — 1 ) • T 2 ! (k — l)(k — 2) 

_ ( a k ~ *) (*k k + l) (a k - k + 2) (t k - k) (t k - k + l) (bfc - k + 2) 

3! (k — 1) (k — 2) (k — 3) 


T J + 


(_ 1 )k-l (&k - *) ( a k — k + l) ... (a k - g) (tt - k) (bk - k + l) . . . (~bfc — g) ^ l 

(k — 1) l (k — 1) ! 
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TECHNICAL NOTE NO . 1716 

TABLES OF HYFERGEOMETRIC FUNCTIONS FOR USE IN 
7 COMPRESSIBLE-FLOW THEORY 
By Yera Huckel 

SUMMARY 


In the holograph method of treating plane potential compressible 
flows the differential equation, originally obtained by Chaplygin in 
his study on gas jets, plays a significant role. This paper tabulates 
various hypergeometric functions which arise as particular solutions -of 
Chaplygin's differential equation. The tables should prove useful in 
the tabulation of other auxiliary functions which may arise in various 
compressible-flow problems. The adiabatic index for air has been 
taken as 1.1*. 


INTRODUCTION 


Any general theory of compressible potential flow will probably 
involve the holograph variables. A reason for this statement is that 
in the holograph plane, in which the independent variables are the 
magnitude and the direction of the fluid velocity, the equations of 
motion are linear; whereas, in the physical plane they are in general 
an intractable set of nonlinear partial differential equations. 

The simplification due to the use of holograph variables, however, 
presents certain difficulties which do not appear in the physical plane. 
For example, application of the necessary boundary conditions for 
uniform compressible flow past an arbitrary body is almost impossible, 
at least up to the present time. Certain singularities in the flow also 
appear, notably, near the sonic speed and in the undisturbed flow at 
Infinity. Nevertheless, the possibility of getting around these and 
other difficulties in the near future Justifies the publication of tables 
for the fundamental set of functions which represent the particular 
solutions of the flow equations in the holograph plane. 

The following section contains equations and definitions necessary 
for the understanding of the several functions listed in the tables. 

The reader is referred to the original paper (reference l) in which the 
particular flow' solutions are derived in detail . 
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EQUATIONS AND DEFINITIONS 


The linear equations in the holograph variables 9 and q, which 
relate the velocity potential (ft and the stream function 1 { f for the 
steady two-dimensional flow of a nonviscous compressible fluid, are 





in which, for the adiabatic equation of state between the pressure and 
density. 


Mi) = 7?! 


a 

(1 - t) P 


and 


Mi> = : iA Ci) 


cLq. v>qy 




1 - (2P + l)r 
q(l - t)P +1 


where 

q magnitude of fluid velocity 

9 angle included Ty velocity vector and positive direction of x-axis 
p density of fluid 
a velocity of sound in fluid 

M Mach number (q/a) 
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7 - 1 


7 ratio of specific heats at constant pressure and at constant volume, 
taken as 1.4 for air 


t dimensionless speed variable ( t = 


2&a. 0 d 20 


M 2 A 

+ M2y 


The index o refers to stagnation point q = 0. 

Observe that the Mach number is given in terms of r by the 
relation 


M 2 = 


1 - T 


For the tables the numerical value 0 = 2-5, corresponding to 7 
is used. Hence 


Mr - 


1 - T 


and M = 1 corresponds to T = g* 

By substituting in equations (l) the product-type solutions 


0* = fkfl) sta («) 


*k = «fc(4) (-**) 


and by observing that from equations (l) 

, n Po clftk(a) 

“k^ - v i-aT - 


a^kti) 


= "kq 


kQt)^ 
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the functions Q^Cq.) can he shown to satisfy the following second- 
order differential equation: 

d2Qt dQv 

q 2 ■—§ + (1 + M2h -f - *2(1 - 112)^ = 0 (4) 

The functions P^Cq.) can he obtained from Qjj-(q) by means of the first 
of equations (3) • Equation (4) may he reduced to a standard type hy the 
introduction -of the dimensionless speed variable t as the independent 
variable. Thus, let 


Qfc(q) = q k Yk( T ) (5) 

where clearly Yjj-(t) 1 as t — ^ 0 (incompressible case) . After 
some elementary operations the desired differential equation is 

t( 1 " T) d^ + [ (k + - ( k + 1 - P) T ] £T + 2 pk(k * l)Yk = 0 (6) 


Equation (6), which is of the hypergeometric type, was first introduced 
hy Chaplygin in his memoir on gas jets (reference 2) . 

In the present- paper, tables of numerical values have been prepared 
for a selected number of the complete set of particular solutions of 
equation (6) . These solutions extend the results of Chaplygin into the 
supersonic range and to negative values of the index k. 


DESCRIPTION OF TABLES 


Tables 1 and 2 have been prepared for the functions Yv and tables 3 

and 4 for the functions r— * for both positive and negative values of the 

dT 

index k ranging from 0-5 to 15 in increments of 0-5 and for the speed 
variable t ranging from 0.01 to 0-50 in increments of 0-01. The 
critical value of t is l/6 for air; hence the present tables extend 
considerably into the supersonic range. Thus, corresponding to the 
value t = 0.50, the Mach number is /~5* 

For large values of* the index k (for example, greater than 15 ) , 
it is possible to develop and utilize asymptotic expressions which 
involve the function 1l(t) for M <1 (see discussion following 
equation (42) of reference l) and involve the function 0(M) for M >1 
(see equation (57) of reference 1) . 
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The numerical evaluation of the functions listed in the tables 
was performed with both manual computing and with the aid of an 
IBM computing machine. The tables may be considered accurate as 
listed although the actual computations made full use of the capacity 
of the machine and involved many more places. 

It is hoped that the tables presented in this paper will be found 
adequate and useful for the numerical evaluation of auxiliary functions 
which may arise in the solution of problems of compressible flow. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 20, 1948 
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